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ABSTRACT 

Considerations of the collision losses for protons traversing the 2.7 K black body 
microwave radiation field have led to the conclusion that the highest energy cosmic 
rays, those observed at > 10 20 eV, must come from sources within the present epoch. 
In light of this constraint, it is here suggested that these particles may be accelerated 
near the event horizons of spinning supermassive black holes associated with presently 
inactive quasar remnants. The required emf is generated by the black hole induced 
rotation of externally supplied magnetic field lines threading the horizon. Producing 
the observed flux of the highest energy cosmic rays would constitute a negligible drain 
on the black hole dynamo. Observations with upcoming air shower arrays and space 
missions may lead to the identification of candidate dormant galaxies which harbor 
such black holes. Although the highest energy events observed so far are accounted for 
within the context of this scenario, a spectral upper bound at ~ 10 21 eV is expected 
since the acceleration to higher energies appears to be precluded, on general grounds. 
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1 INTRODUCTION 

The highest energy cosmic ray induced atmospheric air 
shower yet observed (Bird et al. 1995) corresponds to initia- 
tion by a 3.2 x 10 20 eV particle. Yet, due to energy degrad- 
ing "GZK" interactions with the pervasive microwave back- 
ground (Greisen 1966; Zatsepen & Kuz'min 1966), such a 
baryonic particle {e.g., proton) that energetic could not have 
originated from a distance beyond 50 Mpc; for sources be- 
yond 100 Mpc all proton energies are less than 10 20 eV upon 
arrival here, independent of their initial energy. As pointed 
out by Elbert & Sommers (1995), accelerating such particles 
by powerful radio sources, the presumed best candidates for 
achieving these energies, is not possible for this particular 
event since there are no appropriate radio objects in the cor- 
rect general direction that are within an acceptable distance. 
What other population, one adequately represented locally, 
should be considered as an alternate candidate? Although 
quasars are not a local phenomenon, the present epoch is in 
fact the preferred one for dead quasars (Schmidt 1978). The 
prospect of such apparently inactive quasar remnants serv- 
ing as sources of the highest energy cosmic rays is explored 
within the context of their putative underlying supermassive 
black holes. If the generator driving the high energy particle 
accelerator is a spinning supermassive black hole, the sys- 
tem involved must in all other respects look relatively dor- 
mant {e.g., no obvious jets). If not a dead quasar then what 
other astronomical object could possibly satisfy this require- 
ment? Considering the apparently plausible (albeit specula- 
tive) conditions here identified, this "bottom-up" scenario 



appears to provide a unique basis for a potential explana- 
tion that is remarkably well matched to the puzzle. The 
more speculative highly interesting "top-down" possibility 
disregarded in this present discussion would have the pri- 
mary particles produced at ultra-high energies in the first 
instance, typically by quantum mechanical decay of some 
supermassive elementary particles related to grand unified 
theories (Sigl et al. 1995; Kuz'min & Tkachev 1998). 

2 HIGH ENERGY ACCELERATORS 

In the present epoch, there is a drastic paucity of quasars 
such as the extremely luminous ones {L <; 10 47 ergs/s) evi- 
dent at large redshifts, those with putative black hole nuclei 
having masses ^ 1O 9 M0. Nevertheless, the expected local 
number of dead quasars associated with the same parent 
population (Schmidt 1978; Small & Blandford 1992; Rich- 
stone et al. 1998) is expected to be relatively large. Studies 
of the X-ray sky indicate a pronounced extragalactic cosmic 
background that arises mainly from accretion powered AGN 
emission at previous epochs (Boldt 1987; Fabian & Bar- 
cons 1992). The present-epoch mass density that has thereby 
been built up in the form of supermassive black holes must 
now be substantially more than represented by currently 
active galactic nuclei (Boldt & Leiter 1995). The black hole 
mass spectrum for the nuclei of active Seyfert galaxies only 
extends up to ~ 5 x 10 7 h~ 2 M Q (Padovani, Burg & Edel- 
son 1990), where h = H /{100 km s" 1 Mpc" 1 ). Consider- 
ing a radiative efficiency of 10% for the accretion powered 
bolometric luminosity of AGNs, Chokshi & Turner (1992) 



Mon. Not. R. Astron. Soc, in press 



Mon. Not. R. Astron. Soc. 000, 000-000 (1999) Printed 1 February 2008 (MN plain TgX macros vl.6) 



have calculated the mass built up over all earlier epochs and 
thereby estimated that the expected local mass density of 
compact galactic nuclei is two orders of magnitude greater 
than that accounted for by Seyfert galaxies. They conclude 
that more than half this local density is contributed by in- 
active quasar remnants which are now black holes of mass 
> 1O 8 /i -2 M0; over 10% of this density is associated with 
black holes of mass > 6 x 10 8 h~ 2 MQ. As emphasized by 
Chokshi & Turner (1992), the local universe is expected to 
be well populated by currently inactive remnants of quasars. 
Based on the mass function described by them we have esti- 
mated that, for h ~ 0.5, the minimum number of black holes 
within 50 Mpc having a mass > 4 x 1O 8 M is ~ 40 and that 
there should be about a dozen or more quasar remnant black 
holes of mass > 1O 9 M0 in this volume. These quasar rem- 
nant expectations are consistent with being lower limits to 
the number of corresponding supermassive black holes in- 
ferred from a recent comprehensive study of massive dark 
objects (MDOs) at the centers of 32 nearby galaxy bulges 
(Magorrian et al. 1998). In this connection we note that the 
number of MDOs within 50 Mpc identified in their sample 
as being more massive than 1O 9 M0 is already 8, comparable 
to the total number of Seyfert 1 AGNs out to that distance. 
This is a lower limit to the total number of such supermas- 
sive objects within this volume since their sample of MDOs 
at the centers of nearby galaxy bulges is incomplete, albeit 
sufficiently large for the correlations sought by them (Magor- 
rian 1998). Using the luminosity function for field galaxies 
(Efstathiou, Ellis & Peterson 1988) to estimate the incom- 
pleteness of their sample suggests that the corrected number 
of present-epoch supermassive MDOs could well be at least 
an order of magnitude greater than that so far observed. 

We use the appellation "dead" for the quasar remnants 
considered here in a more complete sense than that which is 
customary in the conventional description of radio galaxies 
as starved/dead qausars (see, e.g., Rees et al. 1982). The 
latter merely implies that the nuclear luminosity of the ob- 
ject is small, whereas the objects we are referring to are 
inconspicuous as regards jets as well. Although local dead 
quasar remnants are manifestly underluminous, their under- 
lying supermassive black holes are likely to be sufficiently 
spun- up (i.e., after their many "Salpeter" time units of ac- 
cretion history; see Thorne 1974; Rees 1997) to possibly 
serve as high-energy accelerators of individual particles. In 
this scenario (c/. Blandford & Znajek 1977) externally pro- 
duced magnetic field lines threading the event horizon of 
such black holes would, by virtue of the induced rotation, 
generate an effective electromotive force characterized by: 
emf oc cBR, where B is the magnetic field strength and 
R is the effective range over which the concomitant elec- 
tric field is applicable. Scaling to the magnitude for this im- 
pressed B field considered by Macdonald and Thorne (1982) 
and taking R w R g (= GM/c ), the gravitational radius, the 
expected value for the emf is then here estimated as 



emf « 4.4 x W 20 B 4 M 9 volts, 



(1) 



where B 4 = B/(10 4 Gauss) and M 9 = M/(1O 9 M ). As de- 
scribed by Znajek (1978) and Blandford (1979), the massive 
black hole in their model behaves as a battery with an emf 
of up to 10 21 volts, comparable to the above estimate. We 
note that radiative losses for particles accelerated in such 
a dynamo are very much larger for electrons than protons. 



Proton energy losses during the phase of acceleration up to 
~ 10 17 eV are due principally to pair production in collisions 
with ambient photons (Blumenthal 1970). For the low lumi- 
nosity objects considered here [i.e., L < 10~ 4 L(Eddington)], 
however, the associated radiation length is larger than the 
gravitational radius (R g — 1.5 x 10 14 Mg cm ) of the black 
hole, for any plausible emission spectrum from the accre- 
tion disk. Hence, radiative cascades (basic to the Blandford- 
Znajek mechanism), such as could attend the process of elec- 
tron acceleration, would not necessarily constitute a compa- 
rable limiting factor in the present scenario for the accelera- 
tion of the relatively few (favorably disposed) protons that 
need to achieve an energy close to that of exploiting the 
full voltage. The dominant energy loss mechanism for pro- 
tons during the final acceleration phase involves photome- 
son production (Blumenthal 1970; Hill & Schramm 1985). 
Although the total cross-section for this increases somewhat 
with energy, its value is always < 200 microbarns at the 
highest energies of interest (Caso et al. 1998). For the low 
luminosity objects considered here, then, the mean free path 
for such inelastic collisions is larger than R g . Thus proton 
energy losses are not expected to be significant in this sce- 
nario. 

This proton accelerator might be an intrinsic feature of 
an underlying dynamo (non jet producing) somewhat differ- 
ent from that of the Blandford-Znajek model or, alternately, 
a subsidiary aspect of that mechanism. The situation we en- 
visage, therefore, is one in which the accelerator is not op- 
erational in the mode in which quasi-steady conversion of 
the hole's rotational energy into that of luminous radio jets 
is possible, but where acceleration of individual protons to 
the full voltage may occur, perhaps sporadically, at those 
episodes (possibly of brief duration) when the emf is not 
shorted out. 



3 REQUIREMENTS 

The air-shower events observed at > 10 20 eV (Cronin 1997) 
correspond to an equivalent omnidirectional cosmic ray flux 
of about one such energetic particle per [(kilometer) 2 - 
decade]. Assuming this flux arises from a local random dis- 
tribution of dead quasars associated with black hole masses 
^ 1O 9 M0 implies that their average luminosity is < 10 42 
ergs/s in these energetic particles. In contrast to the full 
electromagnetic process discussed by Blandford and Zna- 
jek (1977) for effectively extracting much energy (up to a 
rate <; 10 45 ergs/s for a 1O 9 M0 object) the particle accel- 
erator aspect considered here for sustaining this cosmic ray 
output (<; 10 20 eV/particle) would constitute, by itself, a 
relatively negligible power drain on the rotational kinetic 
energy reservoir provided by the putative underlying super- 
massive (<; 1O 9 M0) canonical Kerr hole. Replenishing the 
particles ejected at high energies (> 10 20 eV) would only 
require a minimal mass input; a luminosity of 10 42 ergs/s 
in such particles (if protons) corresponds to a rest mass loss 
< 1O~ 5 M0 in a Hubble time. Since the energy source in this 
model derives from the black hole itself, no further mass in- 
put is needed; the principal input requirement for the am- 
bient plasma is to generate the magnetic field. The ejection 
of the energetic protons considered in this model would con- 
stitute a leakage of electric charge at a rate amounting to 
less than 0.01% of the total effective current flow associated 
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Table 1. Upper limits for cosmic ray generation by nearby massive (> 1O 9 M0) dark objects™. 
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NGC1399 


5.2 


-1.785 


0.028 


0.14 


3.2 


NGC1600 


11.6 


-2.046 


0.051 


0.13 


6.5 


NGC2300 


2.7 


-2.161 


0.067 


0.30 


3.6 


NGC2832 


11.4 


-1.935 


0.040 


0.11 


5.7 


NGC4168 


1.2 


-2.356 


0.105 


0.57 


3.0 


NGC4278 


1.6 


-1.959 


0.042 


0.31 


2.2 


NGC4291 


1.9 


-1.807 


0.029 


0.24 


2.0 


NGC4472 


2.6 


-2.509 


0.145 


0.46 


5.2 


NGC4486 


3.5 


-2.377 


0.110 


0.34 


5.3 


NGC4649 


3.9 


-2.143 


0.063 


0.25 


4.2 


NGC4874 


20.8 


-2.000 


0.045 


0.09 


8.3 


NGC4889 


26.9 


-1.678 


0.021 


0.06 


6.5 


NGC6166 


28.4 


-1.767 


0.027 


0.06 


7.4 


NGC7768 


9.1 


-1.991 


0.044 


0.14 


5.3 


average 


9.3 




0.06 


0.23 


4.9 



"Based on Table 2 in Magorrian et al. (1998). Here, Mg = M MDO /(10 9 M e ), x = M MDO /M Bulge , rh 
[c 2 (dM / dt)] / Lp^d, M — g a l ax y bulge mass loss rate (iWgj/year) per 1O 12 M0, B4 = I?/(10 4 Gauss), (emf)'20 
(era/) /(10 20 volts) 



with the Blandford-Znajek dynamo. 



4 DISCUSSION 

What are the environmental circumstances of the black hole 
nuclei in dead quasars, and are they conducive to sustain- 
ing the magnetic fields needed for this model? We appear 
to encounter a variety of different situations that can ac- 
count for the low luminosity of such systems. Why are they 
so quiescent? As emphasized by Rees (1997), their environ- 
ment could be almost free of gas, so that very little gets 
accreted; in this case the accretion rate is inefficient, in that 
the cooling is so slow (because of the low densities) that only 
a small fraction of the binding energy gets radiated before 
the gas is swallowed, a regime of advection-dominated ac- 
cretion flow (Narayan 1997). The unusually low luminosities 
associated with the supermassive black holes at the centers 
of nearby bright elliptical galaxies (Loewenstein et al. 1998; 
Fabian & Rees 1995; Fabian & Canizares 1988) have been 
explained in terms of advection-dominated flow (Mahadevan 
1997; Narayan 1997), even when there is ample ambient gas. 
Rather than attempting to address the parameter space for 
this model's potential applicability to all dead quasars or 
study possible alternate models for the microphysics of low- 
luminosity AGNs, we have considered in some detail the phe- 
nomenology of the specific massive dark objects (MDOs) in 
the large sample recently reported by Magorrian et al. (1998) 
and have investigated the quite general constraints they im- 
pose on our proposed scenario. Of the 32 nearby galaxies 
with MDOs considered by Magorrian et al. (1998) 14 have 
compact central masses Mm do > 1O 9 M0. In order to obtain 
absolute upper limits on the emf that could be generated by 
each of these putative supermassive black holes we take that 
the maximum magnetic field possible near the event horizon 



would correspond to pressure equilibrium between it and the 
infalling matter, whereby 



(B 4 ) < 3.7m/M 9 , 



(2) 



where rh is the accretion rate in Eddington units, (Lem 
1.26 x 10 47 M 9 ergs/s), defined as 



m = c 2 (dM/dt)/L l 



(3) 



The galaxy bulge mass (M Bulge), as described by Magorrian 
et al. (1998), constitutes the ultimate reservoir fixing the 
maximum accretion rate possible into the central black hole, 
viz: 



m < 4.6 x 10 /z/a;, 



(4) 



where 11 is the galaxy bulge mass loss rate (Mo/year) per 
10 12 M Q and 



X = M M D0/M B ulge 



(5) 



Typically, fi £ 1 (Renzini & Buzzoni 1986; Mathews 1989). 
The maximum possible magnetic field is obtained from equa- 
tions (2) and (4) as 



B 4 <4.1 x 10" 1 (n/xy'M, 



-1/2 



(6) 



The largest possible emf, that corresponding to the max- 
imum magnetic field, is then obtained from equations (1) 
and (6) as 



(em/) 20 < 0.18( M /x) 1/2 M 9 1/2 , 



(7) 



where (em/)20 = em//(10 20 volts). 

The limits on rh, magnetic field and emf obtained from 
equations (4), (6) & (7) are listed in Table 1 for the 14 most 
massive nearby MDOs studied by Magorrian et al. (1998). 
The average maximum magnetic field is 2.3 x lO 3 /! 1 ' 2 Gauss. 
For /i ~ 1 this is substantially less than the magnetic 
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field strength previously considered by Macdonald & Thorne 
(1982), but comparable to more recent estimates of the 
likely strengths of magnetic fields threading the horizons of 
accretion-disk fed black holes (Ghosh & Abramowicz 1977). 
The average limiting emf (4.9 x lO 20 /^ 1 ' 2 volts) is remarkably 
consistent with what would be required to accelerate the 
highest energy cosmic ray yet observed. What is the prospect 
of actually achieving something close to these limiting val- 
ues for several supermassive MDOs? While a well-defined 
specific accretion model that is clearly compatible with the 
constraints and requirements is still to be identified, we find 
that the possibility of such can not be excluded at this stage, 
on quite general grounds. We anticipate that some version of 
the current standard view of advection-dominated flow (Ma- 
hadevan 1997) could be adequate in this respect for many of 
the MDOs considered here. Exploring the possibly relevant 
parameter space for such models is beyond the scope of this 
discussion. However, the margins associated with satisfying 
the constraints listed in Table 1 could well be suggestive 
of what sort of modification of the theory would be needed 
should the more standard models prove inadequate. Based 
on the 14 MDOs listed in Table 1, it appears that the gen- 
eration of cosmic rays more energetic than ~ 10 21 eV is 
precluded; to a precision commensurate with the statistics 
of this sample, such a spectral cut-off can then be taken as 
a prediction of our proposed scenario. 

5 OUTLOOK 

The Akeno Giant Air Shower Array (AG AS A), currently 
the world's largest, has during 1990-1997 accumulated an 
extensive collection of cosmic ray events above 10 18 ' 5 eV, 
461 of which correspond to more than 10 19 eV (Takeda 
et al. 1998). Although only six of these exceed the criti- 
cally important threshold of 10 20 eV, the overall data-base 
already provides evidence, albeit still statistically limited, 
which suggests that the arrival directions and energies are 
most compatible with a scenario in which sources of ultra- 
high energy protons trace the inhomogeneous distribution 
of luminous matter in the "local" present-epoch universe, 
well within 100 Mpc (Medina-Tanco 1999) . If more data at 
the highest energies confirm this preliminary indication it 
would strongly support models based on cosmic ray proton 
acceleration by present-epoch supermassive black holes and 
preclude the need for invoking cosmologically remote sources 
of extremely energetic exotic hadrons whose special proper- 
ties would allow them to traverse large distances (J> 1000 
Mpc) without their energy falling below 10 20 eV (Chung, 
Farrar & Kolb 1998; Farrar & Biermann 1998). The sample 
of about 100 extraordinary events with energy > 2 x 10 20 
eV expected with the upcoming Pierre Auger array (Cronin 
1997) will come from nearby sources and, if protons, will 
point accurately to the directions of origin {i.e., owing to the 
correspondingly large particle Larmor radius in the weak in- 
tergalactic magnetic field) . Candidate galaxies within the ac- 
ceptable pixels would then be searched for stellar-dynamical 
evidence for central supermassive black hole nuclei (such as 
the MDOs discussed by Kormendy & Richstone 1995, Kor- 
mendy et al. 1997 and Magorrian et al. 1998), here taken to 
be indicative of the dead quasar sources of the highest energy 
cosmic rays. If such a correlation is clearly established, and 
the lack of correlation with strong radio sources persists, 



it would imply that the existence of a black hole dynamo 
is not a sufficient condition for the presence of pronounced 
jets. The OWL (Orbiting Wide-angle Light-collectors) space 
borne NASA mission planned for observing, from above, 
those air showers induced by the highest energy cosmic rays 
is anticipated to have the sensitivity for accumulating an or- 
der of magnitude more such events than expected with the 
Auger array (Streitmatter 1998). 
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